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Abstract-- El creciente interés en generadores eléctricos de 
energías renovables conectados a la red plantea la necesidad de 
dispositivos electrónicos de control en tiempo real, con alta 
velocidad de cálculo y precisión, aprovechamiento sofisticado del 
recurso energético y sin distorsionar la red. Los actuales 
dispositivos de lógica programable de alta velocidad y los 
lenguajes de descripción de material con potencialidad de 
lenguajes de programación, permiten abordar este problema. En 
este trabajo se presenta un software programado en VHDL 
diseñado para una máquina eólica con generador trifásico de tres 
niveles, arquitectura espalda contra espalda con un bus de continua 
entre el frente hacia el molino y el frente hacia la red de 
interconexión. Sus características modulares, parametrizables y de 
bajo costo le permiten adoptarse a otras arquitecturas y fuentes de 
energía. Al residir en forma compacta en un único componente 
electrónico puede evitar la perturbación del ruido electromagnético 
de los grandes generadores 


Palabras clave:- AC generators - Field programmable gate arrays - 
Wind power generation. 


I. INTRODUCCIÓN 


A. Objetivo 


Control automático en tiempo real de la generación de 
corriente trifásica de 50 ciclos inyectable a la red de UTE a 
partir de la energía suministrada por un generador eólico. 


B. Contexto 


Este software fue desarrollado dentro de un proyecto 
específico que no se llevó a cabo (una licitación de UTE 
ganada por una empresa que no tuvo continuidad) de un 
parque eólico de características estructurales que explican el 
tipo de solución encarada: 


e 4 aerogeneradores de eje vertical tipo Darrieus en H, 
sección cuadrada de 40 m, con una potencia nominal de 
500 Kw c/u., con 70% de fabricación nacional. Costo 
total presupuestado: 2.5 millones de dólares. 

e Máquina de velocidad variable con un motor asíncrono 
de corriente alterna en jaula de ardilla sin escobillas. 

e Arquitectura de generador de tres niveles, tres ramas, 
espalda contra espalda con un bus de continua entre el 
frente hacia la máquina eólica y el frente hacia la red de 
interconexión. 

e Modulación por puente de doce transistores IGBT cada 


frente, (cuatro por rama) y diodos en paralelo. 
C. Características del software. 


Es una aplicación diseñada para trabajar en tiempo real de 
generación. Permite gran velocidad de cálculo y precisión, 
con reducida carga computacional. Su programación 
modular hace posible un código compacto para que todo el 
programa resida en un solo micro-componente, evitando 
conexiones vulnerables a la perturbación por ruido 
electromagnético del generador. El diseño incluye 


asignación de parámetros constructivos para ser adaptado a 
arquitecturas diferentes. Para ello se ha recurrido a: 


e Programación en VHDL (Very High Speed Integrated 
Circuit - Hardware Description Language, ANSV/IEEE 


1076-1993), lenguaje descriptivo de diseño de 
circuitos tipo ASIC (circuito integrado para 
aplicaciones específicas). Es un lenguaje de 


programación orientada a eventos con potencia similar 
a los lenguajes informáticos procedurales, pero de 
trabajo en simultáneo. 

e Aritmética de punto flotante de doble palabra con una 
biblioteca propia de funciones matemáticas. 

e Cálculo de senos y cosenos por algoritmo CORDIC 
(rotación de vectores), más rápido que por serie de 
Taylor y más compacto que las tablas de líneas 
trigonométricas. Requiere un mínimo de operaciones 
aritméticas elementales. 

e Componente de silicio FPGA (Field Programed Gate 
Array) de última generación. De bajo costo, que puede 
ser grabado en forma serial en producción industrial o 
semi-industrial, permite alojar todo el software y es 
reprogramable. Optamos por el FPGA Spartan 3AN de 
Xilinx, aunque también pueden usarse componentes 
más económicos según la escala de la máquina. 

e Comunicación para su monitoreo por una puerta USB. 

e Modulación de la corriente trifásica por metodología de 
modulación vectorial por ancho de pulso (PWM). 

e Algoritmo de modulación por coordenadas naturales 
que suministra los vectores para cada rama y sus 
tiempos de trabajo en forma directa, simétrica y 
optimizada, minimizando la carga computacional. 
Incluye la gestión de tiempos muertos para evitar 
colisión entre transistores complementarios. 

e Diseño, programación, simulación y prototipación con 
herramientas estándar de acceso gratuito: La principal 
ha sido ISE Design Suite 10.0, que ofrece 
gratuitamente el fabricante, Xilinx. 

e Dicha herramienta permite grabar directamente el 
componente desde un PC por medio de una puerta 
USB, una metodología universal y sencilla. 

e Toda la información técnica de referencia es accesible 
gratuitamente en Internet. 


II. LO QUE HABÍA QUE HACER 


A. Arquitectura del generador 


Partimos de las necesidades de un generador de inducción 
por corriente alterna (ACIM), preferido por su sencillez y 
durabilidad al no tener escobillas que se desgasten ni imanes 
permanentes que lo encarezcan. El sistema de generación de 
velocidad variable consiste en desacoplar la frecuencia de la 
red y la frecuencia mecánica del rotor mediante una 
estructura de conversión CA/CC/CA “Back-to-Back*: dos 


convertidores en oposición conectados mediante un bus de 
corriente continua, con el generador asíncrono doblemente 
alimentado. Los convertidores multinivel permiten aumentar 
la tensión de trabajo y generar una onda de salida de mejor 
calidad, con disminución de la distorsión armónica y de la 
necesidad de filtrado. 


La máquina asíncrona, trabaja por deslizamiento, diferencia 
relativa entre la velocidad del campo magnético y la 
velocidad del rotor. La velocidad de sincronismo es 
levemente superior a la velocidad del rotor, que al girar 
genera una tensión inducida en el bobinado por efecto 
Faraday o efecto generador: en todo conductor que se 
mueva en un campo magnético se induce una tensión. 


Imaginemos un auto en bajada andando en segunda, la 
gravedad empuja al auto más rápido que el motor. En 
nuestro caso es el movimiento mecánico de la máquina 
eólica contra el movimiento del motor eléctrico. El 
deslizamiento es consecuencia de esa diferencia de 
velocidades. La sobre-aceleración del empuje eólico 
produce un voltaje inducido en el bobinado del motor. 


Como sobre la velocidad del viento no tenemos control 
nuestra técnica será actuar sobre el otro componente, 
introducir en la máquina eléctrica una tensión que la haga 
funcionar un poco más despacio y recolectar el exceso de 
voltaje inducido para alimentar la red. Se trata entonces de 
cómo controlar la velocidad rotórica. 


B. El control vectorial. 


El método tradicional de control por potencial y frecuencia 
es ineficaz cuando se producen picos de corriente por 
cambios bruscos de velocidad y par motor. Una respuesta 
rápida exige controlar el par y el flujo en modalidad de lazo 
cerrado por control vectorial o “orientado al campo”, que se 
realiza directamente, consume menos energía, es más 
eficiente, más barato de operar, y reduce también el costo de 
los componentes de nuestra máquina. Como las corrientes 
del rotor no se pueden medir directamente porque gira en 
una jaula de acero sin conexiones eléctricas, el programa las 
calcula indirectamente a partir de los siguientes datos que 
pueden medirse o se conocen: 


e Las corrientes instantáneas de las fases del estator. 
e La velocidad mecánica del rotor. 

e Las constantes constructivas de motor. 

Por otra pare deberemos también tener acceso a: 

e La velocidad del viento. 

e Las tensiones instantáneas de la red. 


Por resultar especialmente clara reproducimos aquí esta 
explicación conceptual tomada de [1] (traducción mía): 


“La clave para entender cómo funciona el control vectorial 
es formarse una imagen mental del proceso de 
transformación de las coordenadas de referencia... 
imaginarnos su funcionamiento desde el punto de vista del 
estator... se aplica al estator una corriente de entrada 
sinusoidal. Esta señal variante provoca la generación de un 


flujo magnético giratorio. Desde una perspectiva fija, los 
vectores de la corriente del estator y del flujo giratorio se 
presentan como magnitudes de la corriente alterna. Ahora, 
en vez de la perspectiva anterior, imaginemos que podemos 
meternos en el interior del motor. Una vez allí, nos 
imaginamos corriendo junto con el bobinado a la misma 
velocidad que el vector de flujo de rotación... ¡Mirando el 
motor desde esta perspectiva, mientras se mantenga en 
condiciones estacionarias, las corrientes del estator 
aparecen como valores constantes, y el vector de flujo de 
rotación está detenido! ...lo que queremos es controlar las 
corrientes del estator para obtener las corrientes deseadas 
en el rotor (que no se puede medir directamente). Con la 
transformación de coordenadas, las corrientes del estator 
se puede controlar como si fuesen valores de corriente 
continua, usando ciclos de control estándar ”. 


El diagrama de la Fig. 1 ilustra la estructura del control. 
Veamos los pasos del control vectorial indirecto. 


1. Se miden las corrientes del estator en dos fases (la tercera 
se infiere por Kirchhoff) con un conversor analógico-digital. 
También se mide la velocidad física del rotor por un 
mecanismo de “encoder” de señales luminosas A, B más una 
señal Z de ajuste de cero. 


2. Estas corrientes se convierten a un sistema equivalente de 
dos ejes, Ia, IB (Transformada de Clarke). 


la = Ib (1) 


I +21 
If = E (2) 


3. Se hace girar este sistema para alinearse al flujo del rotor 
(Transformada de Park) en un ángulo dado por la última 
iteración del ciclo de control, y esto da las variables Id, Iq 
(directa y cuadratura) en el nuevo sistema, que para 
condiciones de estado estacionario serán constantes. 

Id = Ia.cosO +1 P.send (3) 


Iq = —fa.cos0 +16. send (4) 


4. El ángulo de flujo se determina por las variables Id e Iq, 
la magnetización de la última iteración, la velocidad 
mecánica del rotor (calculada a partir de la medida del 
encoder), las constantes constructivas de la máquina, y el 
ángulo obtenido en la última iteración. Las ecuaciones han 
sido convertidas a una formulación adecuada para el cálculo 
computacional. 


ôl = Id — I M,a (5) 
óIM = dI.QT (6) 
IM = IM, +061M (7) 
ólqg = lq. QT (8) 
Slip = 6lq/1IM (9) 
TAS = vel.Q Ang (10) 
0 = 0a + IAS + Slip (11) 


El subíndice a expresa el valor de la magnitud en la 
iteración anterior, el incremento ð refiere a la variación de la 
magnitud en una iteración. En lo demás: 
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Fig. 1. Diagrama simplificado del control vectorial. Para seguirlo de acuerdo al texto, abajo a la izquierda lectura de intensidades la e Ib y 
señales A, B, y Z del encoder, y los sucesivos pasos hasta el Puente Inversor Trifásico, a la derecha. El círculo punteado indica un pulso de 
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reloj diferente, como se explica en la sección “La arquitectura del tiempo”. 


Ol : Incremento de intensidad en una iteración 

IM : Intensidad de magnetización 

QT : Constante constructiva de tiempo retórico 

SlipT : Deslizamiento eléctrico en un período 

IAS : Incremento de la velocidad eléctrica 

QAng: Constante angular (parámetro constructivo = vel.flujo 
* tiempo de conmutación) 


O. Ángulo de flujo. 


5. Entre los valores de cada magnitud medida y su referencia 
respectiva hay una desviación que se usa como feedback para 
determinar la conducta del sistema, en los controladores 
Proporcional/Integral (PI), que aplicaremos en tres casos: a) 
La velocidad del motor contra la referencia obtenida a partir 
de la velocidad del viento, nos da la referencia del par. b) Id 
contra la referencia de flujo, que es una constante, nos da la 
tensión directa a imponer. c) La referencia del par contra lq 
nos da la tensión cuadrática a imponer. Tenemos así dos 
nuevos componentes de nuestro sistema bidimensional 
estacionario. 

Los controladores PI tienen la forma genérica: 

Out = Kp . err + Ki . ferr dt (12) 
Donde Kp y Ki son los factores de términos de error en ese 
ciclo y error acumulado, que son constantes concretas para 
cada una de las aplicaciones a, b y c. 


6. Una vez obtenidos los dos valores de tensión Vd y Vq en el 
sistema bidimensional estático, debemos hacer el camino 
inverso de las transformadas Anti-Park primero, y Anti-Clarke 
después: 


va = Vd cos- Vq sin (13) 
VB= Vdsinð+ Vq cos (14) Anti-Clarke 
Va = V8 (15) 


Vb A v3 Va (16) 
ve -V2 tv3 Va (17) Anti-Clarke. 


2 
Nótese que Va + Vb + Vc=0 


C. Algoritmo de coordenadas naturales. 


Tenemos ahora un vector “fasor” tensión a imponer, (Va, Vb, 
Vc). Aplicamos el algoritmo PWM para obtener los Estados y 
Tiempos de Trabajo que comandan el puente de transistores. 
Nos inclinamos por un algoritmo en coordenadas naturales 
descrito en [2], como “el primer algoritmo de modulación 
vectorial tridimensional para convertidores... usando 
coordenadas naturales. Esta técnica simplifica en gran 
medida la selección del espacio tridimensional donde se 
encuentra un vector de tensión dado. El algoritmo reduce 
drásticamente la complejidad computacional y los cálculos 
para determinar los vectores activos con sus 
correspondientes tiempos que generan el vector de referencia 
de tensión... sin tener que utilizar cálculos trigonométricos, 
tablas o transformaciones de sistemas de coordenadas que 
incrementan la carga computacional ”. 


El algoritmo está convenientemente explicado en ese 
documento, y en la sección siguiente veremos como lo 
adecuamos a nuestro software, pero veamos aquí sus 
fundamentos conceptuales. 


El vector de referencia se encuentra en un espacio 
tridimensional definido por las tensiones Va, Vb, Vc. Cada 
período de modulación genera una secuencia de cuatro 
vectores de estado (los estados posibles del puente inversor) 
ur(t1), u2(t2), us(t3) y us(t4), donde tı, t2, t3 y t4 son los tiempo 
de encendido de esos vectores. 


La entrada del algoritmo es el vector de referencia 
normalizado, al dividirlo por una constante (en nuestro 


software preferimos multiplicar por su inversa) que depende 
del máximo de amplitud que suministra el bus de continua, y 
llevado a números positivos en un rango de O a 2 (para un 
convertidor de 3 niveles). Luego se trabaja con las tensiones 
escaladas y sus partes enteras para encontrar el sub-espacio en 
que se encuentra el vector de referencia. Los sub-espacios son 
tetraedros definidos por los estados y llenan todo el espacio 
de generación. La búsqueda se hace por comparaciones contra 
la ecuación de distintos planos de corte del espacio vectorial. 


Ubicado el vector, se obtiene la composición lineal de los 
estados que definen ese sub-espacio, por los tiempos de 
trabajo correspondientes. Nos da además la secuencia de 
vectores que minimice el número de conmutaciones de los 
transistores, y la dividimos en dos subperíodos simétricos. 
Veamos a continuación como lo implementamos. 


D. La arquitectura del tiempo 


Toda cuestión de tiempo es siempre una relación entre ciclos, 
cuantas veces cabe un ciclo en otro. 


Nuestro problema: generar una triple onda sinusoidal en 
ciclos de 50 Hz, desfasadas en 120". 


Nuestro recurso: Cada onda será producto de 4 transistores 
que trabajan dos a dos generando una tensión -1,0 o 1 por un 
factor de amplitud, durante un tiempo. Eso nos da una 
poligonal de pulsos cuadriculados, una aproximación 
“pixelada” de nuestra sinusoide como en la Fig. 2. Es la 
modulación por ancho de pulso, PWM por su sigla en inglés. 





Fig. 2. Aproximación a la sinusoide por modulación por ancho de pulso. 


La precisión depende de la cantidad de niveles del 
convertidor y de cuan pequeño sea el “pixel”, tiempo de 
trabajo de los transistores, que se prenden y apagan a la orden 
pero que son dispositivos físicos reales con desgaste y 
recalentamiento y tienen un tiempo mínimo de encendido y 
disipación (en nuestro caso en torno a 300 ns). Debemos 
aproximarnos a la meta dentro de esas restricciones porque no 
podemos darles órdenes que no puedan cumplir. 


Nuestra opción fue dividir el ciclo de la sinusoide en 90 
pasos, una conmutación por paso. Cada conmutación “pinta 
un pixel”. Es decir, tendremos un tiempo de conmutación de 
20 milisegundos (ms) dividido en 90 pasos da tiempo de 
222,2... microsegundos (us). 


20 ms / 90 = 222.2.. us 


O lo que es lo mismo, en frecuencia: 
50 Hz x 90 = 4500 Hz = 4.5 KHz 
El ciclo de conmutación es submúltiplo del ciclo eléctrico. 


Esto da una diferencia angular de 4° entre conmutación y 
conmutación. Nuestro cálculo es en radianes y aritmética 
hexadecimal de punto flotante, pero es el dispositivo 
computacional quien carga con el cálculo y nuestra mente la 
que carga con la representación imaginaria del problema. 
Intuitivamente sabemos que llegamos al máximo de tensión 
entre el paso 22 y el 23 (90%) y volvemos a O en el paso 45 
(1809). 


El punto del ciclo en que deberá encenderse el transistor y 
cuánto deberá permanecer encendido son órdenes numéricas 
al conmutador que salen del cálculo hexadecimal. Dos cifras 
significativas en base 16 requiere dividir nuestro tiempo de 
conmutación en 256 intervalos, pero como dividiremos ese 
tiempo por 2 para llegar a una sucesión de estados que sea 
simétrica serán 512 intervalos. Por lo tanto tenemos 222.222 
us / 512 = 434 nanosegundos (ns), o bien 2.304 MHz (4.5 
KHz x 512). 


Nuestra medida del tiempo viene dada por los pulsos de un 
oscilador de cuarzo, que en nuestro caso es de 50 MHz (20 
ns). Hablamos, por tanto, de CANTIDAD DE PULSOS del 
oscilador. El ÚNICO TIEMPO al que tendremos acceso será 
siempre una cuenta de oscilaciones del cuarzo. Cada 11111 
pulsos del cuarzo comenzamos de nuevo las mediciones y 
cálculos necesarios para una etapa de trabajo que dura 222 
us. Pero ese trabajo es una agenda de órdenes que se imparten 
cada 22 pulsos del cuarzo, 440 ns. 


20 ns x 22 = 440 ns 
20 ns x 11111 =222 us 


Los tiempos están redondeados a enteros, porque solo 
podemos contar pulsos del cuarzo en números enteros. 
Descartamos las diferencias fraccionales porque ajustamos 
nuestra frecuencia a la frecuencia real de UTE y la seguimos 
como una consigna. Detectaremos cada comienzo de ciclo de 
UTE y comenzaremos a contar de nuevo. El paso por O (un 
semi-ciclo) lo vemos midiendo la tensión cada 222 us. Dos 
medidas sucesivas de signo diferente indicarán que ha habido 
un pasaje por 0. De esta forma nuestro sistema de conteo se 
ajusta solo. 


Al mismo tiempo debemos fiscalizar la corriente de UTE para 
detectar posibles fallas que nos obliguen a salir de servicio, 
como alteraciones en la frecuencia O la amplitud. Es un 
problema interesante medir una magnitud con una vara de 
medida que cambia, y fiscalizar al mismo tiempo esos 
cambios. La paradoja se resuelve sencillamente matando dos 
pájaros de un tiro, controlando la forma sinusoidal, asumiendo 
que cualquier falla importante en la corriente de UTE alterará 
primero esa forma. El método más sencillo es por medio de la 


derivada en un punto, como vemos en la Fig. 3. Dados dos 
períodos de conmutación separados por un período de ciclo 
sinusoidal, controlamos que el valor absoluto de la relación 
entre las diferencias de tensión en cada período se mantenga 
por debajo de un límite de alarma que a su vez puede 
ajustarse. De esta manera podemos detectar la anomalía antes 
de que ocurra, ganando en seguridad y costos del material 
eléctrico. 
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Fig. 3. Detección de un pico en ciernes por el método de la derivada en un 
punto. (Excluye los períodos en que la derivada es O o cambia de signo). 


Para concluir hablemos brevemente del problema de los 
“tiempos muertos”. Como dijimos, los transistores trabajan de 
a pares: mientras uno esté encendido, su complementario 
deberá estar apagado y viceversa., cambiando de estado 
cuatro veces en 222 us, según órdenes que indican la cantidad 
de pulsos de reloj de 20 ns en que el transistor permanecerá 
en ese estado. 


Es una relación ESTIMULO/RESPUESTA. Pero esto siempre 
implica el problema material de la posibilidad real de la 
respuesta. Si diésemos al mismo tiempo la orden de 
encendido a uno y de apagado al otro, como el cumplimiento 
no puede ser instantáneo ni debemos pretender que lo sea, 
estaríamos generando un estado de colisión. Para evitarlo, 
daremos primero orden de apagado a un transistor y 
demoraremos la orden de encendido a su complementario el 
tiempo suficiente como para que la disipación de la corriente 
haya ocurrido. Esa demora es siempre una cantidad entera de 
intervalos de 20 ns. 


Los transistores son fabricados en serie, pero tienen entre sí 
pequeñas diferencias que determinan su comportamiento, el 
tiempo real de encendido y apagado. Lo encontramos en 
forma empírica en un trabajo de sintonía que se hace en el 
tiempo de montaje de dada máquina eólica concreta. 
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Fig. 4. Demora de encendido del transistor T2 luego de la orden de apagado 
de Tl, en que las curvas representan el comportamiento real de los 
transistores. 


Sacaremos en ello ventaja de la condición reprogramable del 
FPGA. Antes de dejar en él nuestro software de gestión, 
trabajaremos en tiempo de montaje con el FPGA y el puente 
de transistores con un programa específico de monitoreo que 
nos permitirá realizar esa imprescindible tarea: el ajuste. 


IHI. DESARROLLO DEL SOFTWARE. 


A. La herencia semántica. 


Es una ley de la evolución tecnológica que cada innovación se 
exprese en el lenguaje de la tecnología que desplaza. La 
columna habla como árbol, el automóvil como carro de 
caballos, la potencia de un motor de vapor se mide en 
“caballos de fuerza”. 


El FPGA es un dispositivo semiconductor que contiene 
bloques de lógica cuya interconexión y funcionalidad se 
puede programar. Surge como evolución de los dispositivos 
de lógica programada PLD (Programable Logic Device), que 
a su vez desplazan a la tecnología de lógica cableada. Unos y 
otros contienen elementos lógicos programables, pero hay una 
gran diferencia de densidad. Medida en puertas lógicas 
equivalentes un CPLD tiene decenas de miles de puertas, y un 
FPGA cientos de miles, o millones. Además hay una 
superioridad de arquitectura que dota a este último de una 
gran flexibilidad. Un FPGA equivale a muchos CPLD 
conectados y trabajando en conjunto. 


Pero en materia de representación del diseño, la forma en 
como se presenta sigue siendo en gran parte la de la lógica 
cableada. Representamos lo que ocurre dentro del FPGA 
como si fuesen componentes electrónicos separados, unidos 
entre sí por “cables” que no existen. Esto es un recurso del 
diseño, e importa aclararlo para una interpretación correcta de 
lo que sigue. En el esquemático de nuestro software se 
representan como siendo una misma cosa las conexiones entre 
componentes físicos (por ejemplo el reloj de cristal de cuarzo) 
y componentes lógicos que son funciones programadas dentro 
del FPGA, como si fuesen componentes físicos de la 
generación anterior. 


Ahora ya estamos en condiciones de presentar nuestro 
software en concreto, y lo haremos en tres puntos: 

° Organización modular 

° Estructura del código de programación 

e Aritmética y biblioteca matemática 


B. La organización modular 


Nuestro software cuenta con 25 módulos. Presentar todos 
ellos y sus conexiones sería demasiado confuso. Preferimos 
mostrar 4 y explicar la lógica usada. En la Fig. 5 mostramos el 
módulo Control, que comanda la secuencia, los módulos 
Clarke y Park, a modo de ejemplo, que realizan el cálculo de 
las transformadas descriptas, y el módulo Espera que al final 
de la secuencia genera un loop hasta que se cumpla el tiempo 
222 us en que debe iniciarse una nueva secuencia. 


|" Control 








Fig. 5. Esquemático, mostrando 4 módulos y las señales de conexión. 


El módulo Control recibe la entrada directa del reloj de 
cuarzo (clk), y lleva la cuenta de los pulsos. 


Tenemos un bus bidireccional (bidir) en el que viajan todos 
los datos entre módulo y módulo. Así, Clarke pondrá allí las 
variables que resultan de su cálculo, que serán las que toma 
Park para su propio cálculo, y pondrá allí sus resultados para 
el módulo siguiente. 


En el bus bidireccional bidir viajan también las direcciones 
indicando al módulo Control qué módulo de trabajo debe 
entrar en ejecución en ese momento, lo que es devuelto a los 
módulos de trabajo en el bus oe de 8 bits. La señal rw 
(read/write) es activada por Control e indica si el módulo que 
se despierta deberá entrar en modo de lectura (lee los datos 
del bus de datos) o de grabación (los graba allí). Cada módulo 
indica la finalización de cada etapa a Control por medio de 
una señal de salida propia fnn, como lo indican las muchas 
“patitas” de Control que aparecen aquí desconectadas (porque 
los módulos correspondientes no están) y las dos conectadas 
de Clarke (f04) y Park (f05). 


La señal lado es necesaria porque un solo FPGA y un solo 
programa comanda simultáneamente AMBOS puentes 
inversores, el que mira a la máquina eólica y el que mira a la 
red: terminados los cálculos de modulación de un lado 
comienzan los del otro, y se comunican todas las señales para 
los 24 transistores al módulo Conmutador. Pero entre ambas 
secuencias hay algunas diferencias que deben ser tomadas en 
cuenta en cada caso, y el módulo correspondiente toma en 
cuenta ese parámetro en su procedimiento. 


Una vez terminados ambos lados, el programa queda en loop 
entre Control y Espera, hasta que el tiempo del ciclo de 
conmutación se cumpla. Mientras tanto, el módulo 
Conmutador seguirá realizando los cambios de estado en los 
tiempos previstos. 





El bus bidir de 104 bits destina 96 bits para portar 3 datos de 
32 bits (lo que se justifica en la sección de aritmética) y 8 bits 
a direccionar los módulos de trabajo (podemos direccionar así 
hasta 255 módulos). En cambio las señales de fin de tarea de 
cada módulo de trabajo son particulares; las unimos dentro 
del módulo de Control por medio de una sentencia declarativa 
VHDL: 


fre<= f01 and f02 and f03 and f04 and f05 and f06 
and f07 and f08 and f09 and f0a and f0b and f71 
and ff1 and ff2 and ff3 and +ff4 and ff5 and ff6 
and ff7 and ff8 and ff9 and ffe and fff ; 


C. El código de programación 


Cada módulo de trabajo se despierta con el evento de un 
cambio en la señal rw si al mismo tiempo oe señala a ese 
módulo. 


Con rw en O (lectura) almacena los datos del bus bidir en un 
almacenamiento interno (inp), y pone en O (terminé de leer) la 
señal f que le corresponde. 


Con rw en 1 (grabación) realiza los cálculos que le 
corresponden y luego junta la identificación del módulo que 
sigue (sig) y los resultados (out) en el bus bidir, y pone en 1 
(terminé de grabar) la señal f que le corresponde. 


El proceso está esquematizado en el diagrama de la Fig. 6. 


rw=1 y oe=yo0 


rw=0 y Oe=yo0 


CÁLCULOS ->out 


Ff<-0 siggout ->bidir 
f<- 1 


A 


imp <- bidir 


Fig. 6. oe: bus de identificación del mód. de trabajo; yo: cte. que identifica al 
módulo; rw: señal de lectura o grabación; bidir: bus bidireccional de 
direcciones y datos; imp: datos que entran; out: datos que salen; sig: módulo 
que sigue en secuencia. 


El módulo Control ve el estado de los módulos de trabajo en 
el producto lógico fre de todas sus señales de fin. Si fre pasa 
de 0 a 1, un módulo de trabajo acaba de terminar y se deben 
poner los datos en el bus bidireccional y despertar el módulo 
siguiente en lectura; cuando fre pasa de 1 a O quiere decir que 
el módulo completó la lectura, se lo despierta en modo de 
grabación y se libera el bus bidireccional para que el módulo 
pueda grabar. 


Este esquema simple se aplica en los distintos módulos 
difiriendo solamente, salvo algunas excepciones, en la rutina 
concreta de CALCULOS de cada uno de ellos. Las 


excepciones son los casos en que se agregan señales de 
entrada externas u otras señales de reloj, son: Conversor 
analógico/digital, Lector de las señales del encoder, 
Comunicación Ethernet, Conmutador, y Tiempos muertos. 


D. Aritmética y biblioteca de funciones matemáticas. 


El software trabaja en aritmética de punto flotante en doble 
palabra (4 bytes), teniendo 1 byte para exponente y 3 bytes 
para la mantisa: 24 bits, 22% - 1. Los 8 bites del byte de la 
izquierda se destinan: bit 7 al signo del número, bit 6 al signo 
del exponente, y bits 5 a O al exponente. Los números se 
normalizan llevando la primera cifra faccionaria significativa 
a la primera posición luego del punto decimal. 


EXPONENTE MANTISA 


32 24 0 
PD 


Esta precisión es más que holgada para nuestras necesidades, 
pero hemos definido nuestros objetivos de desarrollo más allá 
de la aplicación misma, tratando de ir diseñando herramientas 
de software utilizables en una escala más amplia. 


Biblioteca de funciones matemáticas: El software incluye su 
propia biblioteca, con las siguientes funciones: 


Flotación: Convertir un número a punto flotante. 
Alinea números a igual mantisa para sumarlos. 

Suma algebraica en punto flotante. 

Multiplicación en punto flotante. 

División en punto flotante. 

Rebaja de punto flotante a punto fijo. 

CORDIC: seno y coseno de un ángulo. 

Administrar pila de datos para memoria y recuperación. 


Como dijimos, se buscó definir una biblioteca de funciones de 
uso más amplio, y además permitir la máxima compactación y 
modularización del código. Las funciones son definidas para 
ser invocadas por los módulos de trabajo como Calls en su 
proceso de cálculo usando la misma lógica de 
funcionamiento. Los datos de imput y los resultados se pasan 
por el mismo bus bidireccional de uso general. De esta forma 
el proceso va del módulo de trabajo a la función de biblioteca 
y luego vuelve. Entra varias veces en mismo módulo. 
Usaremos una variable (ope) para identificar la operación en 


curso, cada vez que se cumpla un paso de cálculo haremos: 
ope<=0ope+1; 


Veamos el ejemplo del módulo que realiza la transformada de 
Park que implica varios pasos de cálculo. Omitimos el código 
de cada rama del case. 


case ope is 


when 1 => -- Recuerda FI anterior 

hen 2 => -- Rebajar FI anterior 

Te Es -- Cordic: seno y coseno de FI 
batas -- IAlfa * Coseno 


when 5 => -- IBeta * Seno 

Meis -- Alinea p. sumar en flotante 
TRN -- IAlfaxCoseno + IBetaxSeno 
es -- -IAlfaxCoseno + IBetaxSeno 
A => -- terminé el cálculo, 


---sigo con próximo paso 


end Case; 


Algoritmo CORDIC. Un caso particular de nuestra 
biblioteca matemática es el módulo que realiza el cálculo de 
seno y coseno de un ángulo (ambas líneas simultáneamente) 
por el algoritmo CORDIC de rotación de vectores 
(Coordínate Rotation Digital Computer), que se usa 
ampliamente con este tipo de recursos computacionales 
porque las únicas Operaciones son sumas, restas y 
corrimientos binarios (multiplicaciones y divisiones por 
potencias en base 2). El procedimiento consiste en obtener las 
coordenadas x e y en cada aproximación a partir de la 
anterior. 


x = X + yAnterior * signo(z) * 2.; 
y = y - xAnterior * signo(z) * 2.; 
Zz = Z - signo(z) * arctan(2.i) 


El algoritmo wusa los arcotangentes de las potencias 
fraccionarias de 2, o sea arctg( 1⁄2) , de 1⁄4, etc.). Definimos 
esas constantes en una tabla dentro del código VHDL de 
nuestra función: 


constant atg00 
constant atg01 
constant atg02 
sawet 


: integer := 16#c910# ; 
: integer := 16#76b2# ; 
: integer := 16#3eb7# ; 


-- arctg(1) * 65536 
-- arctg(1/2) * 65536 
-- arctg(1/4) * 65536 


IV. RESULTADOS EXPERIMENTALES 


El desarrollo del software se hizo en forma completa en fase 
de programación, puesta a punto y simulación en las 
siguientes etapas: 


° Desarrollo en MATLAB y SIMULINK del algoritmo de 
modulación PWM, y del algoritmo CORDIC, y prueba. 

e Desarrollo de cada módulo VHDL y prueba por 
separado en OrCAD y en ISE 10.0, por medio de un 
Test Bench para cada uno. Lo mismo en las distintas 
secciones multi-módulo del software, como ser: lectura 
del encoder, transformadas de los sistemas de 
coordenadas, pruebas de las funciones matemáticas, etc. 
para ir integrando las partes progresivamente. 

e Integración completa del software y Test Bench integral 
en ISE 10.0 en Behavorial Simulation. 

° Síntesis y compilación, y Test Bench funcional y de 
timing. 

a Por último la etapa de Fit (grabación física) del software 
en el FPGA y prueba en funcionamiento real del mismo 
testeando su salida desde la placa por medio de un 
osciloscopio, comprobando la señal sinusoidal de salida. 
Para ello se desarrolló también un módulo de auto- 


testing integrado y grabado junto con el software 
aplicativo. 


Lo que seguía luego era la prueba de la placa en un prototipo 
de generador eólico en escala. Esa etapa no pudo ser 
cumplida por el cierre de las actividades de la empresa. 


IV. PERSPECTIVAS 


Este es un software para un proyecto de mega-generación de 
un parque eólico de unas pocas torres grandes. En ese 
contexto era más práctico y económico instalar directamente 
el Kit de Desarrollo de XILINX para el Spartan XC3S700AN, 
con un costo de USS 200 FOB. Esta solución puede adoptarse 
para cualquier proyecto de mega-generación, sin problemas. 


En otros contextos será necesario analizar en concreto el 
diseño de la placa más conveniente a implementar y la forma 
de grabación en serie. Eso depende del costo integral de la 
placa en relación al costo final de la máquina eólica y de la 
escala de producción, según se trate de mega, meso o micro- 
generación. El estudio de factibilidad deberá ser hecho en 
cada contexto. 


Un relevamiento primario de los pasos subsiguientes para la 
implementación del proyecto nos lleva al siguiente marco: 


e El componente FPGA Spartan XC35700AN tiene un costo 
FOB por separado de U$S 37, y representaría el costo 
principal de una placa integral que tuviese además: 

o memoria NOR FLASH paralela de 4 Mb que carga 
el programa 

puerta USB 

reloj de cuarzo de 50 Mh 

dos conversores A/D con un reloj independiente 

los conectores FX2 de I/O para los drivers de 

transistores y otras señales 

microcontrolador para manejar la USB 

o puerta RS232 


O O O O 


O 


e MONTAJE y TESTING: Debe resolverse el tema del 
encapsulado y la inspección de soldadura. Aunque puede 
hacerse localmente en pequeñas cantidades, ya para el 
orden de 1000 unidades se hace necesario orientarse a la 
fabricación en el exterior. 


e Como alternativa el software puede alojarse también en 
componentes más baratos, caso del XC38200A. (U$S 13) 
y con un encapsulado bastante más amigable. 


El proyecto debe estudiarse concretamente en la escala de la 
máquina que se quiera montar. La conveniencia está 
obviamente en la diferencia de costos entre la importación de 
componentes con montaje nacional y fabricación de las partes 
en que haya capacidad de hacerlo, contra la importación del 
puente inversor ya completo. En esa ecuación de costos, el 
costo del software que se está pagando junto con el 
inversor importado hace una gran diferencia. 


V. CONCLUSIONES 


Presentamos aquí un software desarrollado en el marco de un 
proyecto concreto de mega-generación eólica. El software 
está orientado a satisfacer las necesidades de control en 
tiempo real, con alta velocidad de cálculo y precisión, y 
aprovechamiento inteligente del recurso energético para 
lograr una mejor relación costo/beneficio y la calidad 
necesaria de la onda suministrada de modo de no distorsionar 
la red de interconexión eléctrica con perturbaciones. 


Para estos objetivos se buscó incorporar tecnología avanzada 
tanto en el desarrollo de software como en la integración de 
los componentes electrónicos en los que va grabado para la 
finalidad específica. El estudio se orientó a esa necesidad, 
optándose por la programación completa en VHDL y la 
grabación en un componente FPGA reprogramable de última 
generación y bajo costo. El proyecto incluyó el desarrollo de 
una biblioteca propia de funciones matemáticas y otras 
herramientas aplicables en una escala más amplia. 


El objetivo de desarrollo de software fue cumplido a 
satisfacción y en un plazo total (incluyendo la etapa de 
investigación) de algo menos de año y medio. 


Aunque el proyecto marco no pudo continuar por razones de 
otra índole, la estrategia de desarrollo adoptada permitió 
llegar a un producto versátil que puede ser aplicado, tal cual o 
con adaptaciones, en una escala mucho más amplia dentro de 
este tipo de aplicaciones hoy en pleno desarrollo tecnológico 
y comercial. 


Como conclusión, podemos decir que la posibilidad de 
incorporación de tecnología de punta o los costos de su 
desarrollo NO SON EL OBSTÁCULO en esta área 
productiva. 
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